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INTRODUCTION 


Automated  structural  optimization  methods  using  finite-element  analysis 
are  now  available  with  sufficient  capability  to  be  used  in  practical  design 
studies.  To  examine  their  use  in  an  actual  design,  the  optimization  of  two 
components  of  a  shuttle  orbiter  vehicle  is  described.  The  computer  program 
used  is  the  Automated  Structural  Optimization  Program  (ASOP)  described  in 
Reference  1.  Several  extensions  of  this  program  which  were  developed  to 
increase  its  generality  are  also  described. 

Many  different  orbiter  and  booster  configurations  were  examined,  during 
Grumman' s  initial  studies  of  the  space  shuttle  system.  In  these  studies  the 
ASOP  program  was  used  extensively  on  the  wing  and  tail  surfaces.  However, 
because  of  the  lack  of  experience  with  the  optim'  _,n  of  fuselage  structures, 
no  automated  optimization  of  the  fuselage  was  initially  attempted.  For  the 
study  reported  herein,  use  of  the  program  for  optimization  of  both  the  wing 
and  the  fuselage  was  undertaken  to  demonstrate  the  versatility  of  ASOP. 

Results  obtained  from  the  automated  design  procedure  are  compared  with  the 
structural  design  obtained  by  more  traditional  methods. 

The  particular  design  chosen  for  this  study  is  known  as  the  H-3T 
configuration;  its  general  structural  arrangement  is  shown  in  Figure  1.  The 
orbiter  vehicle  for  this  design  has  a  delta  wing  and  two  main  propulsion 
engines.  Four  air  breathing  engines  are  also  provided  for  flight  within  the 
atmosphere.  The  liquid  hydrogen  for  the  main  propulsion  engines  is  carried 
in  two  external  disposable  tanks  attached  to  the  side  of  the  fuselage.  The 
crew  compartment  is  also  jettisonable  in  an  emergency  and  is  not  structurally 
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part  of  the  fuselage.  Access  to  the  payload  is  through  clam  shell  doors 
along  the  top  of  the  fuselage.  These  doors  are  also  not  part  of  the  primary 
structure. 

This  report  shows  how  both  the  wing  and  fuselage  structures  are  idealized 
into  finite  element  mathematical  models.  This  will  serve  to  guide  the  program 
user  in  establishing  the  idealizations  needed  for  the  optimization  process. 

To  demonstrate  the  program's  usefulness,  results  obtained  from  the  automated 
calculations  for  the  element  gages  are  compared  to  material  distr:'  '■‘•ions 
obtained  by  the  more  traditional  non-automated  methods.  To  demonstrate  the 
program’s  convergence  character'  sties,  the  fuselage  structure  was  optimized 
twice  starting  with  different  material  distributions,  with  essentially  the 


same  final  results 


Figure  1  General  Structural  Arrangement  of 


DISCUSSION 


A.  Wing  Design  Study 

1.  Structural  Idealization 

The  general  planform  of  the  wing  is  shown  in  Figure  1.  The  semi- 
span  from  the  centerline  of  the  fuselage  to  the  tip  is  about  650  inches; 
the  root  chord  is  approximately  800  inches.  The  development  of  the  finite- 
element  model  of  the  primary  structural  box  evolves  naturally  from  the  defined 
positions  of  the  spars  and  ribs  as  shown  in  Figure  2.  Each  of  the  covers, 
outboard  of  the  root ,  is  idealized  as  anistropic  membrane  quadrilateral 
and  triangular  elements  which  carry  all  shear  and  axial  stress  -  i.e.,  no 
bar  elements  are  used  ( see  Appendix  B) .  The  webs  of  the  spars  and  ribs  are 
treated  as  shear  panels  separated  by  posts  which  enable  transfer  of  external 
loads  to  the  structure.  Inboard  of  the  root,  the  type  of  construction  dictates 
the  need  for  introducing  bar  elements  to  represent  the  carry-through  beams. 
Sufficient  attention  to  the  details  of  the  carry-thi-ough  structure  is  required  to 

properly  account  for  the  boundary  effects  on  the  outer  portion  of  the  wing. 

The  configuration  employs  covers  of  stringer-reinforced  titanium 
sheet.  The  internal  construction  uses  standard  titanium  spars  and  ribs.  3y 
using  the  spar  and  rib  inter  section  as  the  nodal  points  for  the  finite  element 
model,  the  exposed  wing  is  broken  into  ten  spanw>se  segments  and  an  average 
of  five  segments  in  the  chvrdwise  direction.  Since  the  distribution  of 
pressure  in  the  chordwise  direction  is  expected  to  be  smooth  and  since  there 
are  no  internal  cutouts  in  the  covers,  the  subdivisions  in  the  chordwise 
direction  should  be  adequate  to  represent  the  true  state  of  stress.  All  the 
elements  in  the  cover  of  the  wing  are  well  shaped,  that  is,  the  triangular 


UPPER  COVER 


Orbiter  Wing 


elements  are  approximately  45°  isosceles  triangles  and  the  quadrilateral 
elements  are  all  of  aspect  ratio  of  two  or  less  with  the  exception  of  the 
last  element  in  the  trailing  edge  root  region.  This  element  is  not  connected 
to  a  fixed  point  on  the  boundary  and  the  strain  energy  in  the  element  can  be 
expected  to  be  low.  For  this  reason  the  discrepancy  between  the  predicted 
and  true  stress  distributions  can  be  expected  to  be  unimportant  to  the  overall 
behavior  of  the  structure. 

The  spars  and  ribs  are  represented  by  shear  panels.  The 
axial  load-carrying  capacity  of  the  web  material  near  the  covers 
is  assumed  to  be  accounted  for  in  the  adjacent  cover  elements,  a 
factor  which  must  be  considered  in  the  final  interpretation  of  the 
results  retained  in  the  computer  program. 

In  constructing  the  geometry  data  for  structures  where  the  primary 
stresses  arise  from  bending  moments,  for  example  wings  and  tails,  a  pre¬ 
liminary  estimate  must  be  made  of  the  radius  of  gyration  so  that  the  behavior 
of  the  structure  will  be  accurately  reflected  in  the  finite  element  model. 
Consider  the  portion  of  a  wing  cross  section  shown  in  Figure  3. 


£  Centroidal  Axis 


Figure  3  Portion  of  a  Typical  Wing  Cross  Section 
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Here  t  is  the  spacing  between  stringers,  b  is  the  depth  of  the  stringer  from 
the  center  of  the  cover  sheet,  and  x  is  the  height  of  the  cover  sheet  above 
the  centroidal  axis  of  the  wing.  The  moment  of  inertia  of  the  finite  element 
model  about  this  axis  should  be  the  same  as  the  actual  wing,  that  is: 

I  =  V  A  =  ♦  £  a/  ♦  a/  ♦  A>t(x  -  §) 

where  y  is  the  height  of  the  element  above  the  neutral  axis  of  the  finite  element 
model,  A  is  the  total  cross  section  area  of  the  cover  sheet  of  width  l  and 
thickness  t.  I  is  the  moment  of  inertia  of  the  stringer,  A  is  the  sheet  area 

XX  3 

between  stringers  and  A  ^  is  the  area  of  the  stringer  (it  is  assumed  that  the 
centroid  of  the  stringer  is  at  b/2  from  the  center  of  the  cover  skin).  Making 
the  cross  section  area  of  the  model  and  the  actual  structure  the  same  and  defining 
R  as  the  ratio  of  stringer  area  to  sheet  area,  one  obtains 

x  =y2A.IK.+  i  As  +  As  x2  +  E  (X  -  f)' As 


2  2  12 

or  neglecting  b  compared  to  x  and  thereby  neglecting  1^,  and  ^  Agt  as  compared 

to  the  remaining  terms  in  the  expression  for  inertia, 
y2  (1  +  R)  =  x2  (1  +  R)  -R  b  x 

2  2  Rbx 

°r  y  =  x  -  jjfej  t 


Here,  y  is  the  desired  height  at  which  the  finite  element  nodal  point  should  be 
above  the  centerline.  The  thickness  of  the  elements  is  now  t*  =  t  +  the 

average  thickness  of  the  cover  and  the  element  idealization  would  be  as  shown 
in  Figure  4. 
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Figure  4  Finite  Element  Idealization  of  Typical  Cover  Element 


We  have  now  replaced  the  actual  structure  by  a  finite  element  idealization 
that  will  have  the  same  bending  stiffness.  One  should  note  that  there  is  a 
slight  decrease  in  torsional  stiffness  because  the  two  covers  are  now  closer  to¬ 
gether,  and  there  is  a  slight  increase  in  torsional  stiffness  because  the  covers 
contain  material  that  is  normally  in  the  stringers.  In  general,  if  these  effects 
are  regarded  as  significant,  the  covers  should  be  represented  by  a  combination 
of  bars  and  membrane  elements. 

2.  Loading  Conditi.riS 

Because  of  the  preliminary  design  nature  of  these  studies,  exact 
pressure  distributions  on  the  wing  were  not  available.  Consequently, 
simplified  but  realistic  distributions  had  to  be  assumed.  Five  loading 
conditions  were  selected  as  being  critical.  These  include:  two  maximum 
dynamic  pressure  conditions  at  positive  and  negative  angles  of  attack;  a  2,5g 
pullout  condition  on  descent  from  orbit; and  two  landing  conditions.  The 
magnitudes  of  the  loads  are  shown  in  Figure  5.  The  distribution  of  the  loads 
on  the  exposed  surface  of  the  wing  was  assumed  to  be  as  shown  in  Figure  6. 

These  distributions  consist  of  a  subsonic  configuration  and  a  transonic 
configuration.  The  subsonic  distribution  peaks  rather  sharply  at  the 

leading  edge  and  tapers  to  zero  at  the  trailing  edge,  with  the 
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Spanwise  Distribution: 


Subsonic 


Chordwise  Distribution: 


Subsonic  ra 


Transonic 


Figure  6  Distribution  of  Applied  Loads  on  the  Orbiter  Wing 
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spanvd.se  distribution  decreasing  uniformly  from  root  to  tip.  The  transonic 
distribution,  in  contrast,  has  a  linear  distribution  frcm  leading  to 
trailing  edge  and  a  constant  magnitude  from  root  to  tip.  These  distributions 
are  adjusted  to  give  the  correct  total  load  and  center  of  pressure  locations 
for  each  of  the  loading  conditions.  The  determination  of  the  loads  applied 
to  a  particular  nodal  point  is  based  on  the  pressure  acting  at  that  nodal 
point  and  its  associated  cover  area. 


3.  Design  Criteria 

The  material  assumed  for  the  optimization  is  titanium  Ti-6Al-6V-2Sn 
annealed  having  the  room  temperature  properties  given  below: 


II 

fi 

155 

ksi 

<< 

ii 

145 

ksi 

F  = 

148 

ksi 

cy 

F  = 

100 

ksi 

su 

E  = 

17.0  X  10  xsi 

In  using  the  ASOP  program,  curves  of  allowable  stress  versus  the  appropriate 
structural  index  must  first  be  developed,  tabulated,  and  then  applied 
to  account  for  the  instability  of  the  cover  panels  (see  Reference  2).  The 
allowable  stress  curve  used  for  this  particular  example  is  the  curve  for 
Y-stiffened  panels  shown  in  Figure  7.  A  length  of  thirty  inches  is  assumed 
between  supports. 
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Conventional  Cover  Allowable  Compressive  Stress 


4.  Optimization  Results 

The  wing  idealization  was  optimized  using  the  fully  stressed  design 
option  of  the  ASOP  program.  Three  iterations  were  performed.  After  three 
iterations  the  element  sizes  change  only  in  the  third  significant  figure, 
Therefore,for  this  idealization  no  more  iterations  are  necessary.  The  top 
cover  thicknesses  are  shown  in  Figure  8.  The  thicknesses  increase,  as  expected 
from  tip  to  root.  The  only  exceptions  are  at  a  few  locations  near  the  tip 
where  the  local  applied  loads  are  high  and  the  wing  is  still  fairly  thin, 
resulting  in  high  bending  stresses.  The  thicknesses  of  the  spars  and  ribs 
are  almost  all  at  the  minimum  of  .02  inches  and  for  this  reason  are  not  shown. 
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B.  Fuselage  Design  Study 


1.  Structural  Idealization 

The  idealization  of  the  orbiter  begins  at  fuselage  reference  station 
4o8,  (See  Figure  1  ).  Forward  of  this  point  the  structure  is  essentially 
an  aerodynamic  fairing  and.  is  not  considered  primary  structure.  The  bulk¬ 
head  station  408  is  treated  as  a  simple  combination  of  bars,  beams,  and  shear 
panels  as  shown  in  Figure  9  .  The  shear  panels  and  beams  allow  the  shear, 
axial  load  and  moment  from  the  forward  tip  of  the  fuselage  to  be  introduced 
into  the  structure  as  concentrated  loads  at  this  station.  The  bulkhead 
and  frame  idealizations  at  each  of  the  important  fuselage  stations  is  shown 
in  Figures  10  and  11.  It  should  be  understood  that  the  idealization  includes 
lumped  properties  of  -he  adjacent  members  in  the  real  structure  and  not  only 
the  properties  of  the  particular  bulkhead  or  frame  alone.  Some  specific 
comments  on  these  structural  components  are  now  given. 

The  next  bulkhead  in  the  idealization  is  at  station  548.  The 
corresponding  structure  is  shown  in  the  structural  arrangement  drawing  in 
Figure  11.  The  important  point  to  note  when  comparing  the  idealization 

and  the  actual  structural  arrangement  (Figure  11  )  at  this  station  is  the 
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level  of  detail  preserved  in  the  idealization.  The  number  of  shear  panels 
used  in  the  bulkhead  is  the  minimum  necessary  to  preserve  the  size  and 
location  of  the  cutout  for  the  forward  oxygen  tank  support,  and  to  provide 
the  backup  structure  for  the  external  hydrogen  tank  support.  There  are 
two  beam  elements  to  maintain  the  approximate  external  contour  across  the 
top  of  the  bulkhead . 

Another  point  of  interest  is  station  608.  This  is  the  point  of  attach¬ 
ment  of  the  cabin  to  the  fuselage.  The  bar  element  representing  this 

attachment  can  be  seen  in  the  idealization.  The  frame  is  composed  only  of 
beams  around  the  outside  and  bar  elements  on  the  interior.  The  bar  elements 
are  needed  to  support  node  52  and  to  represent  the  lumped  areas  of  the  shear 
panels  between  the  adjacent  bulkheads. 

The  next  bulkhead  at  station  686  serves  mainly  as  a  point  of  application 
for  a  portion  of  the  distributed  loads  on  the  fuselage  and  also  serves  to 
break  up  the  shear  panels  along  the  side  of  the  fuselage  into  elements  of 
reasonable  aspect  ratio. 

Station  825  is  at  the  bulkhead  between  the  cabin  and  the  cargo  bay. 

Because  the  cargo  bay  covers  are  non-structural,  the  bulkhead  is  idealized 
only  up  to  the  hinge  line  of  the  cover.  The  aft  cabin  support  is  at  node  83. 
This  dictates  the  form  of  the  subdivision  of  the  bulkhead  into  shear  panels. 
Again  the  outside  of  the  bulkhead  is  defined  by  beam  elements. 

The  bulkheads  at  stations  955  and  1083  are  similar.  Each  is  a  series 
of  beam  elements  around  the  fuselage  countour  and  bar  elements  defining  the 
central  keel  of  the  fuselage  and  the  floor  of  the  cargo  bay.  In  this  region 
the  fuselage  is  essentially  a  two  cell  box  beam.  The  structural  idealization  at 
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these  two  stations  can  be  compared  to  the  typical  construction  in  this  region 
shown  in  the  structural  arrangement  drawing  at  station  1138  (Figure  11). 

The  next  section  of  interest  is  at  station  1213.  This  bulkhead  supports 
the  aft  end  of  the  main  oxygen  tank  and  receives  the  most  forward  wing  spar. 

This  is  also  the  beginning  of  the  air  breathing  engine  bay.  The  bar  from 
node  159  to  173  is  added  to  stabilise  one  of  the  fuselage  longerons  that  ends 
at  this  bulkhead  and  it  is  not  structural.  At  this  station  the  bulkhead  has 
again  been  divided  into  as  few  shear  panels  as  is  necessary  to  preserve  the 
basic  behavior  of  the  bulkhead. 

The  next  idealized  bulkhead  is  at  1326  where  the  air  breathing  engines 
are  connected  to  the  structure.  The  bulkhead  at  station  ibOO  is  the  aft  booster 
attachment  bulkhead.  At  this  point  the  booster  thrust  loads  are  transferred 
to  the  orbiter  making  the  longeron  just  forward  of  the  interstage  attachment 
point  one  of  the  heaviest  members  in  the  structure. 

The  bulkhead  at  station  1502  provides  aft  support  for  the  external 
hydrogen  tank  attachment.  The  beam  elements  that  make  up  the  perimeter  of 
the  frame  portion  do  not  represent  actual  members  in  the  structure,  but  are 
the  result  of  lumping  together  the  properties  of  the  frames  on  both  sides  of 
station  1502. 

The  thrust  structure  for  the  main  orbiter  propulsion  engines  is  provided 
by  bulkheads  at  1572  and  1622.  The  structural  arrangement  drawing  (Figure  11 ) 
is  presented  only  at  1572  which  is  also  the  end  of  the  cargo  bay  and  the  forward 
attachment  point  for  the  vertical  fin. 

There  are  several  interesting  details  in  the  structural  idealization  of 

this  bulkhead  that  should  be  noted.  One  of  the  shear  panels  in  the  original 
idealization  had  to  be  replaced  by  a  bar  element  because  the  adjacent  sides  of 


18 


the  element  were  almost  parallel.  The  main  thrust  structure  is  represented  by 
a  horizontal  beam  between  1572  and  1622.  The  forward  side  of  this  beam  is 
represented  by  the  bars  and  shear  panels  between  nodes  273  and  257,  and 
274  and  259.  The  shear  webs  for  the  beam  are  represented  by  shear  panels 
between  the  two  bulkheads.  Because  this  bulkhead  represents  the  end  of  the 
cargo  bay, support  must  be  provided  for  the  shear  panels  comprising  the  inner 
shell  of  the  bay.  These  supporting  nodes  and  the  shape  of  the  inside  of  the 
cargo  bay  can  be  seen  by  tracing  out  the  nodes  274,  283,  28l  and  279. 

The  bulkhead  idealization  at  1622  shows  the  rear  end  of  the  thrust  struc¬ 
ture  beam.  Also  attached  at  this  bulkhead  is  the  support  for  the  on  orbit  LOg 
tank  which  occupies  most  of  the  space  above  the  center  of  the  bulkhead. 

The  last  two  bulkheads  are  almost  rings  because  of  the  large  cutouts  for 
the  rocket  exhaust.  The  tops  of  these  two  are  attachment  points  for  the  fin. 
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Figure  10  continued.  Fuselage  Bulkhead  Idealizations 


sr/ix&a 


ST/U8Z5 


Figure  11  Fuselage  Bulkhead  Structural  Arrangements 
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Figure  U  continued  Fuselage  Bulkhead  Structural  Arrangements 
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Figure  11  continued  Fuselage  Bulkhead  Structural  Arrangements 
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Figure  11  continued  Fuselage  Bulkhead  Structural  Arrangements 
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Figure  11  continued  Fuselage  Bulkhead  Structural  Arrangements 
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2.  Loading  Conditions 

Sixteen  loading  conditions  were  selected  as  possibly  being  critical  for 
different  portions  of  the  orbiter  fuselage.  A  schematic  of  the  mission 
profile  is  shown  in  Figure  12.  The  points  where  the  different  loading 
conditions  occur  are  indicated.  A  summary  of  these  loading  conditions,  with 
the  magnitude  of  the  safety  factor  for  each  condition,  is  given  in  Figure  13# 

The  magnitude  of  the  loads  in  each  of  the  conditions  and  other  relevant 

infornation  is  given  in  Figure  l4  .  To  make  this  information  more  meaningful, 

the  envelopes  of  ultimate  bending  moment,  shear,  and  axial  loads  are  plotted 

in  Figures  15  ,  16,  and  17  .  In  these  figures  the  critical  loading  condition 

for  each  point  along  the  fuselage  is  also  indicated.  As  can  be  seen,  the  critical 

bending  moment  conditions  over  much  of  the  fuselage  is  load  condition  5  during  boosted 

flight  when  the  booster  thrust  line  is  several  feet  below  the  orbiter  centerline. 

The  interface  loads  between  the  wing  and  fuselage  and  between  the  vertical 
fin  and  fuselage  were  obtained  by  calculating  the  reaction  loads  of  the  wing  and 
tail  when  optimized  independently  as  cantilevered  structures.  For  a  more 
rigorous  analysis  of  these  interactions  it  would  be  necessary  to  do  a  coupling 
analysis  of  the  wing,tail  and  fuselage  or  alternatively  analyze  the  entire 
vehicle  as  one  structure.  The  interaction  loads  for  the  orbiter  and  booster 
were  obtained  from  a  previous  analysis  of  the  two  vehicles  together. 
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Figure  12  Mission  Profile 
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Figure  14  Summary  of  Applied  Loading  Conditions 


Orbiter  Fuselage  Ultimate  Bending  Moment  Envelope 


3.  Design  Criteria 


The  orbiter  fuselage  is  assumed  to  be  made  of  the  same  material  as  the 
wing  as  given  previously.  But  because  in  the  fuselage  large  portions  of  the 
skin  are  in  diagonal  tension  the  elastic  properties  of  most  of  the  shear 
panels  must  be  modified.  The  only  exceptions  are  several  panels  near  the  main 
propulsion  engines  which  are  heavily  loaded  and  therefore  thick  enough  to  resist 
shear  buckling. 

To  modify  the  shear  stiffness  of  the  panels  in  diagonal  tension,  the 
shear  modulus  was  reduced  to  75 %  of  its  normal  value.  The  justification 
for  this  is  based  on  the  graph  in  Figure  18  which  is  copied  from  Reference  3* 


Figure  18  Modulus  Ratio  for  Buckled  Elastic  Web 
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In  this  figure  t/t  *s  the  ratio  of  shear  stress  to  shear  buckling  stress. 

Gjdj/G  is  the  ratio  of  effective  shear  modulus  for  the  buckled  panel  to  shear 

modulus  for  an  unbuckled  panel.  A^^/dt  is  the  effective  area  of  the  upright 

(in  this  particular  case  a  ring)  to  the  cross  section  area  of  the  sheet.  For 

a  wide  range  of  t/t  ,  the  shear  stiffness  ratio  is  about  .75  when  the  upright 
•  cr 

area  is  large  compared  to  the  sheet  area  between  uprights.  In  the  ASOP 
program  the  input  is  in  terms  of  E  (Young’s  modulus)  and  v  Poisson's  ratio 
and  it  is  therefore  not  possible  to  specify  G  directly.  Instead  orthotropic 
material  properties  are  specified  for  all  shear  panels  in  diagonal  tension. 

The  failure  criteria  used  for  the  fuselage  are  the  standard  criteria  in 
ASOP,  that  is, maximum  shear  stress  in  the  shear  panels  and  maximum  allowable 
axial  stress  in  the  bar  elements  that  are  buckle-resistant.  The  buckling 
stress  criteria  used  for  optimization  for  those  bar  elements  that  are  not  well 
supported,  for  example  the  hydrogen  tank  attachment  struts, is  shown  in  Figure  19- 
4.  Optimization  Results 

Two  separate  optimizations  of  the  fuselage  were  performed  starting  with 
different  initial  element  sizes.  The  initial  optimization  started  with  a 
size  distribution  that  was  generated  by  hand.  First,  the  structure  was  sized 
using  simple  strength  of  materials  formulas  and  the  ultimate  applied  moment 
and  shear  envelopes.  Using  these  rough  sizes  a  finite  element  analysis  was 
performed  to  obtain  accurate  internal  loads.  With  these  internal  loads  the 
structure  was  again  resized  by  hand.  The  weight  of  the  resulting  finite 
element  model  was  17,200  poundr.  Except  for  a  few  minor  exceptions,  all  the 
shear  panels  in  the  structure  were  at  the  minimum  allowable  thickness.  The 
structure  was  then  optimized  using  ASOP.  Three  iterations  were  performed  using 
the  fully  stressed  design  option  in  the  program.  The  weights  of  the  structure 
for  the  iterations  were  12,150  lbs., 12,200  lbs.  and  12,080  lbs.  The  major 
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difference  between  the  final  structure  weighing  12,080  lbs.  and  the  initial 
structure  weighing  17,200  lbs.  was  the  areas  of  the  frame  elements  which  were  not 
initially  sized  for  the  applied  loads,  instead  only  approximate  values  were 
used.  The  most  interesting  results  of  the  optimization  are  the  sizes  of  the 
main  longerons.  For  convenience  these  longerons  have  been  identified  in 
Figure  20-  The  areas  of  these  longerons  at  the  start  of  the  optimization 
and  after  the  final  iterations  are  given  in  Figure  21.  Also  the  critical 
loading  condition  is  indicated.  As  can  be  expected  from  an  examination  of 
ultimate  bending  moment  envelope,  load  condition  5  is  the  critical  condition 
over  most  of  the  length  for  those  longerons  near  the  top  and  bottom  of  the 
fuselage . 

In  general,  the  longerons  near  the  forward  and  aft  ends  of  the  fuselage 
are  at  minimum  area.  From  a  comparison  of  the  starting  size  r;nd  the  final  size 
it  can  be  seen  that  the  initial  calculations  using  strength  of  materials 
theory  and  a  single  finite  element  analysis  gave  a  material  distribution  that 
is  fairly  close  to  the  final  distribution.  This  is  important  because  it 
indicates  that  the  ASOP  program  is  generating  material  distributions  similar 
to  those  that  would  be  generated  by  more  traditional  methods. 

After  this  first  optimization  of  the  structure  was  completed,  another 
one  was  performed  starting  with  uniform  element  gages.  The  results  of  this 
exercise  are  shown  in  Figure  22  where  the  longeron  areas  are  given  in  a  form 
similar  to  Figure  21.  This  time  the  areas  for  the  first  resizing  and  the 
last  resizing  are  given.  Seven  iterations  were  performed  in  all.  The 
critical  loading  conditions  are  the  same  as  before  and  are  not  repeated 
in  the  figure.  Hie  weight  of  the  structure  after  three  iterations  is  essentially 
the  same  as  the  final  weight  in  the  previous  case.  The  weight  after  seven 
iterations  is  11,750  pounds.  An  examination  of  these  final  element  sizes  and 
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Figure  22  Areas  for  the  H-3T  Longerons  After  Optimization 
Starting  with  Uniform  Sizes 


the  final  sizes  for  the  previous  example  shows  that  the  two  designs  are 
essentially  the  same  even  though  the  starting  points  were  different. 
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APPENDIX  A 


Modifications  for  Anisotropic  Load  Distribution 


Equation  4-2  pages  48  and  49  of  reference  1  has  been  modified  to 
give  the  correct  load  distribution  for  anisotropic  construction.  This 
has  been  accomplished  by  using  the  values  of  E  in  the  direction  of  the 
element  edges  in  formula  4.2.6. 


This  is  done  by  pre  and  post  multiplying  the  elastic  constant  AiJ  by  the 


transformation  matrix 


2  2 
cos  y  |  sin  y 


"*■  “2 


,  2 

-sin  v 
sin  2y  >  -sin  2y 


sm  y  »  cos  y 

cos  y  i  sin  Y  *  cos  Y 

I  “ 

cos  2y 


where  y 


is  the  angle  between  the  elastic  property  axis  and  the  side 
Tliis  gives  a  new  elastic  property  matrix  A  .  Prom  A^,  E 


of  the  element. 

can  be  easily  extracted. 
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APPENDIX  B 


ORTHOTROPIC  BEHAVIOR  OF  INTEGRALLY  STIFFENED  PANELS 


A  non-isotropic  structure  frequently  encountered  in  present  day  aircraft 
structural  analysis  is  the  integrally  stiffened  panel.  Fig.  23.  The  following 
derivation  is  made  with  the  purpose  of  applying  the  anistropic  finite  elements 
herein  described  to  problems  of  analysis  of  such  structures,  where  the 
stiffeners  form  a  regularly  distributed  orthotropic  gridwork.  This  work  does 
not  consider  the  flexural  behavior  of  such  structures  but  does  describe  the 
in-plane  mechanics  of  deforuiation. 


b  =  Stiffener  spacing  in  x  direction 
a  =  Stiffener  spacing  in  y  direction 

=  lijt^  area  of  stiffener  in  x  direction 

Aj>  =  h^tg  area  of  stiffener  in  y  direction 

A 

q  =  -L  _  Stiffener  area  -  x  direction 
B1  at  ~  Plate  area  -  x  direction 


-?  .  -  1  +  fifCl^y) 

JL 

V  x  + 


Stiffener  area  -  y  direction 


Plate  area  -  y  direction 


Consider  the  strain  energy  of  the  module  pictured  on  P.  47 
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Where  oj  ,  {^  ,  xKy  ,  £x  ,  £y  ,  Yty  are  stress  and  strain 
components  in  the  plate.  G]  ,  6 1  are  stress  and  strain  of  the  stiffener  in 
the  x  direction  while  ^  ,  £2  are  the  stress  and  strain  of  the  stiffener 
in  the  direction. 


Assumptions ; 


1.  Stress,  and  therefore  strain,  does  not  vary  with  coordinates 
(i.e.  independent  of  x,  y). 

2.  For  compatibility  of  the  structural  components  we  must 
assume  fc,  =  €x  and  £  t  =  £  y 


Stress-Strain  Relations: 
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Substituting  the  stress-strain  relations  into  the 
strain  energy  equation  and  performing  the  indicated  volume  integrations, 
we  obtain  the  following  expression  for  strain  energy. 


In  order  to  obtain  generalized  stress  components  we  must 
define  a  strain  energy  density  function  with  respect  to  the  plate 
volume  in  order  that  integration  over  the  plate  volume  yields  the  total 
strain  energy  in  the  stiffened  panel. 
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It  must  be  noted  that  the  modified  Hooke’s  Law  that  follows 
is  based  on  this  step.  When  the  modified  Hooke*s  Law  is  used  to  derive 
a  stiffness  or  flexibility  matrix,  the  integration  must  be  done  over 
the  plate  volume  only.  The  strain  energy  density  is  therefore: 


()  -  '  *1  £* 


»  +  cr  *  ~ +G>  <T* 


The  generalized  stress-strain  relations  are  obtained  by 
differentiating  the  strain  energy  density  with  respect  to  the  components 
of  the  strain  tensor,  for  example; 


€  = 
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The  resulting  modified  stress-strain  relations  are: 
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From  assumption  (2)  and  equations  (b)  and  (h)  we  can  obtain  the 
actual  stress  in  the  plate  and  ribs. 
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